NEUTRON STAR ENVELOPES AND THERMAL 
RADIATION FROM THE MAGNETIC SURFACE 



JOSEPH VENTURA 

Physics Department, University of Crete, and lESL, FORTH 
71003 Heraklion, Crete, Greece 

AND 

ALEXANDER Y. POTEKHIN 

loffe Physico-Technical Institute, 194021 St. Petersburg, Russia 



Abstract. The thermal structure of neutron star envelopes is discussed with 
emphasis on analytic results. Recent progress on the effect of chemical constitution 
and high magnetic fields on the opacities and the thermal structure is further 
reviewed in view of the application to pulsar cooling and magnetars. 



1. Introduction 

Neutron stars (NS) are formed with very high internal temperatures ap- 
proaching 10^^ K in the core of a supernova explosion (see, e.g., Shapiro &: 
Teukolski 1983 - ST83 in the following). Copious neutrino emission brings 
the temperature in the stellar core down to ~ 10^ K within about one day 
and then, more gradually, to ~ 10® K within 10"^ years. It was realized early 
on that such objects were likely to have effective surface temperatures of the 
order of 10^ K (Chiu & Salpeter 1964). Comparison with theoretical cooling 
curves can further provide information on aspects of the internal structure 
of NS such as superfluidity and the possible appearance of pion or kaon 
condensates and strange matter in their interior (e.g., Pethick 1992; Page 
1997; Tsuruta 1998). It is then clear that the observation of this surface 
cooling X-ray emission is an objective of prime scientific importance. 

The observation of such faint point sources has turned out to be diffi- 
cult however, having to await the modern era of imaging X-ray telescopes. 
ROSAT observations in the 90's finally yielded improved spectral informa- 
tion, opening a new chapter in our ability to probe the internal structure 
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of superdense matter. Joachim Triimper (this meeting) already gave us an 
overview of the recent observations in this field (see also Caraveo et al. 
1996; Triimper & Becker 1997). Increasing interest is also presently be- 
ing directed on the new and important subject of magnetars, extensively 
discussed during this meeting. 

Theoretical models of NS cooling are of necessity rather complicated, 
requiring detailed understanding of the stellar structure, equation of state 
(EOS), and thermal balance over an enormous range of densities and chemi- 
cal diversity (e.g., Pethick 1992; Tsuruta 1998, and references therein). Dur- 
ing an initial period of 10^ — 10^ y the star cools principally via neutrino 
emission from its interior. There are many neutrino emission mechanisms, 
whose rates depend on the state of matter, particularly on nucleon super- 
fluidity (for a review see Yakovlev et al. 1999). An older NS cools mainly via 
photon emission from its surface. Finally, the cooling may also be influenced 
by heating processes such as friction due to differential rotation between the 
superfluid and normal parts of the star (e.g., Alpar et al. 1989), /3-processes 
arising from chemical imbalance during the spindown (Reisenegger 1997), 
and pulsar polar cap heating due to impinging charged particles accelerated 
in the magnetosphere (e.g., Halpern & Ruderman 1993). 

Analysis of the NS thermal evolution is considerably simplified however 
by the fact that the stellar interior, from densities of ~ 10^° g cm~^ inward, 
is nearly isothermal because of the very high thermal conductivity in these 
layers. It is therefore convenient to establish a relation between this interior 
temperature Tj, defined as the temperature at mass density p = 10^° g 
cm~^, and the effective surface temperature (Gudmundsson et al. 1983 - 
GPE in the following). Thus, it becomes possible to examine separately the 
properties and structure of the outer envelope, which turns out to be crucial 
in determining the ratio Tg/Ti and the nature of the emitted radiation (e.g., 
GPE; Hernquist & Applegate 1984 - HA84 in the following; Ventura 1989; 
Potekhin et al. 1997 - PCY in the following). 

In the next section, we consider basic features of the mechanical and 
thermal structure of the outer NS envelope without magnetic field. The 
strong magnetic field (B ^ 10^^ G), which was found in most of the pul- 
sars, shifts the atmosphere bottom and the region of partial ionization to 
higher densities. Furthermore, it strongly affects the radiative and thermal 
transport through the envelope. These effects will be discussed in Sect. 3. 

2. Outer Envelope of a Neutron Star 

The outer envelope consists mainly of electrons and ions. It extends down 
to a depth of a few hundred meters, where the density 4.3 x 10^^ g cm~^ 
is reached. At this density, neutrons begin to drip from the nuclei (e.g.. 
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ST83). Thus the inner envelope, which extends deeper down to the core, 
consists of electrons, atomic nuclei, and free neutrons. 

The outer envelope can be divided into the atmosphere, the liquid ocean, 
and the solid crust. The outermost layer constitutes a thin (0.1 - 100 cm) 
NS atmosphere (optical depth r < 1), where the outgoing radiation is 
formed. The plasma density at the atmosphere bottom is about 0.001 - 0.1 
g cm~^, depending on temperature, surface gravity, and chemical compo- 
sition. This plasma can be partially ionized and non-ideal. Bound species 
can be distinct until the electron Fermi energy becomes comparable with 
the Thomas-Fermi energy at p < 10 Z A g cm~^, A and Z being the mass 
and charge numbers. At higher densities the ions are immersed in a jellium 
of degenerate electrons, which still strongly responds to the Coulomb fields 
of the ions as long as p < lOZ'^A g cm""^. 

In the rest of the envelope (at p^ 10 Z'^A g cm~^) the electrons form a 
strongly degenerate, almost ideal gas. This gas is non-relativistic at pg ^ 1 
and ultrarelativistic at pQ 3> 1, where pe = p/10^ gcm~^. The ions form a 
Coulomb gas or liquid at F < 175, where F ^ 22.75 {pq/ /le)^^^ Z^/^ /Tq is 
the Coulomb coupling parameter, Tq = T/10^ K, pe = A/Z. At T 175, 
the liquid freezes into a Coulomb crystal. The pressure is almost entirely 
determined by degenerate electrons and thus independent of T, while the 
mass density is mostly determined by the ions. 

While cooling of the NS during the initial neutrino dominated era is not 
influenced by the outer layers, it is in fact the properties of these surface 
layers that characterize the flux and photon spectrum emitted at the NS 
surface, leading in turn to estimates of Ti. At the subsequent photon cool- 
ing stage, the heat insulation by the envelope controls the cooling rate. The 
radiation is reprocessed in the atmosphere, which yields an emitted spec- 
trum, in general different from that of a black body (Romani 1987). The 
properties of these layers are thus crucial to interpreting observations and, 
understandably, a lot of theoretical work has been devoted to analyzing 
the thermal structure of non-magnetic and magnetic NS envelopes (e.g., 
GPE; PCY; Heyl &; Hernquist 1998) and the radiation properties of their 
atmospheres (e.g., Pavlov et al. 1995; Page & Sarmiento 1996; Rajagopal 
et al. 1997; Potekhin et al. 1998, and references therein). 

2.1. MECHANICAL STRUCTURE 

To review the cardinal properties of the NS surface layers, let us recall the 
enormous gravitational potential GM/R 0.14&{M /Mq) R^^ (? (where 
Rq = R/10^ cm ^ 1 and M/Mq ~ 1.4 for most typical NS), which renders 
effects of General Relativity appreciable. Indeed, the Schwarzschild radius 
Tg = 2GM/c^ fs 2.95 (M/Mq) km is not much smaller than the stellar 
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radius R (throughout this review, M is the gravitational mass of the star, 
which is 10-15% smaller than its baryon mass). The hydrostatic equilibrium 
is then governed by the Oppenheimer-Volkoff equation (e.g., Thorne 1977). 
Introducing the local proper depth z = {R — r) (1 — rg/R)~^^'^ (where r is 
the radius), in the surface layers {z <^ R) one can rewrite this equation in 
the Newtonian form 

dP/dz = gp, (1) 

where 

g = ^ 1.327 X 10^^ (1 - r,/i?)-i/2 ^ Rf cm s-^ 

is the local gravitational acceleration at the surface. 

Since the surface gravity is huge, the atmosphere's scale height is rather 
small. In the non-degenerate layers, we have 

P = (p//xmu) kT, (2) 

where mu = 1.6605 x 10~^^ g is the atomic mass unit and fi = A/{Z + 
1). Thus, following a thin non-degenerate atmosphere of a scale height 
P/{9P) ^ 0.626 (Tg/^) i?g (M/Mq)"^ cm, electron degeneracy sets in at 
densities /? > 6 /icTg^^ g cm~^. 

The electron kinetic energy at the Fermi surface is kTp = ep — nieC^, 
where k is the Boltzmann constant, 

Tf = (7f - 1) niec^/k = 5.93 x 10^ xV(l + 7f) K (3) 

is the Fermi temperature, and 

jF = \Ji + X, X = = 1.009 — , (4) 

V irieC irieC \fleJ 

are the electron Lorentz factor and relativity parameter, respectively. 

Elementary fitting formulae to the pressure of fully ionized ion-electron 
plasmas as function of density at arbitrary electron degeneracy and temper- 
ature have been presented by Chabrier & Potckhin (1998). These formulae 
can also be extended to partially ionized atmospheric layers in the mean-ion 
approximation, provided the effective ion charge Z is known (cf. PCY). At 
sufficiently high density, however, where Tp ^ T, the main contribution is 
that of strongly degenerate electrons with the pressure depending only on 
the density through the parameter x (e.g., ST83): 
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where Pq = nieC^/X = 1.4218 x lO^^ dynes cm ^ is the relativistic unit of 
pressure, X = ?i/{mec) being the Compton wavelength. This may further be 
approximated as Pg ~ ^'oX^''^'''*/(9vr^7ad), where 7ad is the adiabatic index, 
equal to 5/3 at x 1 and 4/3 at x ^ 1- 

Note that for strongly degenerate electrons 



dP = riedep = 4.93 x 10^^ erg g ^ (p//Lte) xdx/7F- 
Prom Eqs. (1) and (6), one obtains 

m.„r^ 4930 cm 



1.027 ^ = X' 



' z 






(-i)l 







Zo 



nicC 
rriugHe 



l^e 914 



(6) 



(7) 



where gi4 = g/{W cms^ ), and zo is a depth scale at which degenerate 
electrons become relativistic. 

Let us note also that the mass AM contained in a layer from the surface 
to a given depth z is solely determined by the pressure P{z) at the bottom 
of the given layer: 



AM{z) = A-kB? P{z) g-^ {l-Vg/ Rf'^ ^ 1.192x10-^ g^l MP {z)/Po. (8) 



2.2. THERMAL STRUCTURE 

From Eqs. (7) and (8), we see that the outer envelope is a very thin layer - 
typically within the outer 100 m, containing ~ 1O~^M0 , - which renders the 
thermal diffusion problem essentially plane-parallel and one-dimensional. 
Assuming a constant heat flux throughout the outer envelope, the temper- 
ature profile can be obtained by solving the heat diffusion equation: 

dT IQaT^dT 16aT^ 

where F is the heat flux, k is thermal conductivity, a is the Stefan-Boltz- 
mann constant, and K is the usual Rosseland mean over the energy spec- 
trum of the specific opacity. This leads immediately to a temperature pro- 
file T/Tg ~ (|t + ^)^/^, where the local effective surface temperature Tg is 
defined through F = aT^ and the integration constant corresponds to the 
Eddington approximation (r = | at the radiative surface, where T = T^). A 
more accurate boundary condition requires solution of the radiative transfer 
equation in the atmosphere. A distant observer would infer from the spec- 
trum and flux the redshifted surface temperature = Te (1 — rg/RY^"^ 
and apparent radius i?oo = -R(l — rg/R)'^/"^ (e.g., Thorne 1977). 
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A knowledge of the mean opacity K{p, T) is then needed to relate the 
temperature to the other plasma parameters, which can then also be ex- 
pressed as functions of the optical depth r. It is also needed in order to 
compute the overall depth and the temperature ratio Ti/T^. 

2.2.1. Opacity 

Heat is transported through the envelope mainly by radiation and by con- 
duction electrons. In general, the two mechanisms work in parallel, hence 

K = Kr + Kc, K-^ = K-^ + K-^, (10) 

where k^, i^c and Kj-, denote the radiation and conduction compo- 
nents of the conductivity and opacity, respectively. Typically, the radia- 
tive conduction dominates the thermal transport (^r > Kc) in the outer- 
most non-degenerate layers of a NS, whereas electron conduction domi- 
nates (kc > Hr) in deeper, mostly degenerate layers. In the absence of in- 
tense magnetic fields, modern cooling calculations (e.g., PCY) make use of 
the Livermore library of radiative opacities OPAL (Iglesias & Rogers 1996), 
which also provides an EOS for the relevant thermodynamic parameters 
at p < 10 Tg gcm~^. For the electron conduction regime, modern opacities 
have been worked out by Potekhin et al. (1999a). 

Radiative opacities. In order to derive an analytic model of the NS en- 
velope, HA84 and Ventura (1989) have written the atmospheric opacity in 
the form 

K{p,T) = KoP^tP. (11) 

In particular, this relation describes the opacity given by the Kramers for- 
mula, which corresponds to a = 1 and (3 = —3.5. In a fully ionized, non- 
relativistic and non-degenerate plasma, the opacity provided by the free-free 
transitions is (e.g., Cox &; Giuli 1968) 

^ 75 5efr [Z/iil) p Tg-^-S cmVg, (12) 

where p is in g cm~^ and ^eff ~ 1 is an effective dimensionless Gaunt fac- 
tor, a slow function of the plasma parameters. For a colder plasma, where 
bound-free transitions dominate over free-free ones, the Kramers formula 
remains approximately valid, but the opacity K is about two orders of mag- 
nitude higher. An order-of-magnitude (within ~ 0.5 in logi^) approxima- 
tion to the realistic OPAL opacities for hydrogen in the range of parameters 
Tq ~ 10-1 - 10°-5 and p ~ (10^^ - 10^)T^ is given by Eq. (12) if we put 
formally ^eff ~ p"^'^- An analogous order-of-magnitude approximation to 
the OPAL opacities for iron at Tq ~ 10° - lO^'^ and p ^ (10"^ - 10"^) T| 
is obtained with ^eff ~ 70p~°'^. These approximations also belong to the 
class of functions (11), but with a = 0.8. 
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Conductive opacities. Thermal conduction of degenerate matter in deeper 
layers is dominated by electrons which scatter off ions. This conductivity 
can be written as (Yakovlev &; Urpin 1980) 

_ TTfc^rmeC^x' ^ 2 S X 10^5 Te erg 

12Ze4A7| AZTT^^^^TTk' ^^^^ 

where A is the Coulomb logarithm. Accurate analytic fitting formulas to A 
as function of p and T have been obtained recently by Potekhin et al. 
(1999a). In the solid crust, this function is reduced to small values by 
quantum and correlation effects in Coulomb crystals. However, we shall 
see shortly that in not too cold NS, the thermal profile is mainly formed 
in the liquid layers of the envelope. Therefore, for our purpose, it will be 
sufficient to note that in the NS ocean A is a slow function of the plasma 
parameters and can be approximated by a constant of the order of unity. 

In the case of non-degenerate electrons, the conductivity can be found, 
e.g., by the method of Braginskii (1957), which yields 

nf^bx 10^° (Fz/A) Tg/^ erg/(cm s K), (14) 

where Fz is a slow function of Z: for example, i<26 = 1.34 and Fi = 0.36, 
whereas the Coulomb logarithm A is ~ 1 near the onset of degeneracy 
and logarithmically increases with decreasing density. Effectively, Eq. (14) 
may be viewed as an analog to Eq. (13) where the dimensionless Fermi 
momentum x has been replaced by an appropriate thermal average (oc VT) . 

2.2.2. Temperature Profile 

The simple functional form of the opacity allows now an analytic treatment 
of the thermal structure (Urpin & Yakovlev 1980; HA84; Ventura 1989). 

Non- degenerate regime. Using Eqs. (1), (11), and (2) one may rewrite 
Eq. (9) as 

Ap_'^9P.rr._^(i 9 T^dT _16 g ( k yr^+^dT 
F - 3 K Ti "3 KoUmu; ' 

which is readily integrated from the surface inward to give the temperature 
profile. In the region far from the surface, where T ^ T^, the integration 
constants may be dropped. Using again Eqs. (2), one can present the result 
as a simple power law: 

^ 3_ A + a-/3 k KoTl 
po+i iQ a + 1 iivriu g 
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This result depends on our implicit assumption that the effective ion charge 
Z remains constant, which in general is not strictly valid. 

Interestingly enough, as noted by HA84, Eq. (16) establishes that the 
conductivity k is constant throughout the radiative non-degenerate layer: 

16 (tT^ 16 a ^^a-(ikaT^ 

^ — — — [ 1:^ (17) 

Kp 3 Ko /0"+^ 1 + a /U"iu5 ' 

The constant value depends on the emitted flux, but is independent of Kq. 
It follows immediately from Eq. (9) that temperature grows linearly with 
geometrical depth. Furthermore, from Eq. (16) we obtain the T, p profile. 
In particular, substituting a = 1 and = —3.5, we obtain 

K,= }^>^{^)\ „ = 2X10"^^. (18) 

cm Tee \Tc) ^514 cmsK 

Substitution of K from Eq. (12) yields 

Te « 0.284^14 /xz « (SO^eff qf'^'' {p/l^cY'^\ q = T^e ^/(/iffw)- (19) 
In these relations, z and p are measured in CGS units, and Tee = Te/10^ K. 

Radiative surface. It is further interesting to note that the opacity K is 
also slowly varying in this region. Combining Eqs. (11), and (19), we obtain 
K cx p/r3-5 oc p-i/i'^ Invoking Eq. (2), we get K oc p-^/i'^. This justifies 
the assumption Ky- ~ const, which one often employs when determining the 
radiative surface from equation (e.g., GPE; PCY) 

(^rP)surface = (2/3)5. (20) 

Using Eqs. (2) and (12), we obtain 

1 /2 

Ps ~ 0.1 /xe gu) T^t g cm-^'. (21) 

Substituting ^eff — 1 and ^eff — 200 for hydrogen and iron, respectively, 
we obtain ps ~ 0.07 1/514 ^fe'^ gcm~^ and ps ~ O.OOA y^guT^i^ gcm~^ for 
these two elements, in reasonable agreement with numerical results of PCY. 



Onset of degeneracy. The solution given by Eq. (19) can be extended 
down to a depth where the electrons become degenerate. Let us estimate 
this depth from the condition kT^ = pp/2me. We obtain 

Xd - 0.053 {g,s qf'\ Id ^ 8.5 x 10^ {g^s qf'^ K. (22) 

Even for very high ^ 10^ K, we have Xd ^ 1) the electrons are non- 
relativistic at the degeneracy boundary. With decreasing T^, the quantities 
Xd and Td decrease, i.e. the boundary shifts toward the stellar surface. 
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Sensitivity strip. Numerical cooling calculations (GPE) revealed early on 
that accurate knowledge of the opacity law is especially important within a 
certain "sensitivity strip" in the (p, T) plane. The temperature ratio T\/T^ 
changes appreciably if K is modified by, say, a factor 2 within this narrow 
strip, while comparable changes of K outside the strip would leave this 
ratio unaffected within a high degree of accuracy. 

The importance of the layer in the outer envelope where the opacity (10) 
turns from radiative to conduction dominated is now easy to demonstrate. 
In the non-degenerate radiative part the integral over p is dominated by 
the higher densities near the base of the layer, while in the degenerate, 
conductive layer it is dominated by the top, least dense part of the layer. 
The region where K turns from radiative to conductive is thus seen to 
contribute most of the resistance to heat flow. The line in the (p, T) plane 
where Ky- = is easily determined from Eqs. (12) and (13): 



Using again the solution (16), we find explicitly the temperature Tt and rel- 
ativity factor xt at the turning point from radiative to electron conduction: 



Some caution is necessary here, however, because the approximation K ~ 
Ky- is not justified as we approach the turning point: actually, at this point 
K = Kr/2, as seen from Eq. (10). In addition, the extrapolation of the so- 
lution (19) to the turning point is, strictly speaking, not justified, since the 
electron gas becomes degenerate: xt > Xd for most typical NS parameters. 
Nevertheless, since xt and Xd are not very much different, the section of 
the thermal profile where our assumptions are violated is relatively small, 
so Eq. (24) provides a reasonable approximation. This is confirmed by a 
direct comparison with numerical results (PCY), which reveals an error 
within only a few tens percent, provided T > 10^-^ K. 

Solution beyond the turning point. An analytic solution to the thermal 
profile in degenerate layers of a NS envelope has been first obtained by 
Urpin & Yakovlev (1980), based on the conductivity in the form (13). The 
hydrostatic equilibrium of the degenerate surface layers is determined by 
Eqs. (1) and (6), which yield gniviiie = 'TieC^(x/7F) dx/d2;. Using the heat 
diffusion equation (9) and Eq. (13), we obtain 




Tt«i2.3x 10^3f/''(?6/^7 K, 



(24) 



dT _ 12 FZe'^K 7f 
dx TT m-ak'^cgiie X^ 



(1.56 X lO'^ Kf 



ZAT^JF 
l^egu 2x^ 



(25) 
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Treating A, A, and Z as constants, we can integrate this equation from xt 
inwards and obtain 

T\z) = t2 + (1.56 X 10^ Kf [fix) - /(xt)] , (26) 

where /(x) = ln(x + 7f) — Tf/Xj ^^'^ the dependence of x on z is given by 
Eq. (7). 

Let us use the above solution to evaluate Tj. Since xt ^ Ij but x 3> 1, 
it is easy to see that [/(x) - /(xt)] ~ Xt^ in Eq. (26). Using Eq. (24), we 
obtain 

1/2 /7T4\^°/^^ 

Ti=(T^ + Ti]\ rA«4xlOM 26 K, (27) 

where we have neglected some factors which are close to unity. 

Figure 1 illustrates the accuracy of the above analytic solution as well 
as its limitations. Here, solid lines show temperature profiles of a NS with 
mass M = 1.4 Mq and radius i? = 10 km obtained by solving the radia- 
tive transfer equation with spectral OPAL opacities in the atmosphere and 
integrating the thermal structure equation (9) with accurate radiative and 
conductive Rosseland opacities inwards in the deeper layers (Shibanov et 
al. 1998). Dashed curves depict the above analytic approximations. The 
left panel corresponds to an envelope composed of iron, and the right panel 
shows the thermal structure of an accreted envelope with its outermost lay- 
ers composed of hydrogen, which is further burnt into heavier elements (He, 
C, Fe) in deeper and hotter layers. This thermo- and pycnonuclear burning 
is responsible for the complex shape of the upper profile. Different straight 
lines show the points at which thermal profiles corresponding to various 
heat fluxes would cross the radiative surface, the region of the onset of de- 
generacy, turning points Kc = K^, and (on the left panel) the bottom of 
the ocean. The latter line is not present on the right panel, because freezing 
of hydrogen and helium is suppressed by large zero-point ionic vibrations 
(cf. Chabrier 1993). 

One can see that not only the above analytic approximations correctly 
describe the qualitative structure of the envelope, but they also provide a 
reasonable quantitative estimate of the temperature at a given density. The 
accuracy deteriorates for lower Te (especially when Z is high), because in 
this case radiative opacities are affected by bound-bound transitions and 
strong plasma coupling effects and thus no longer obey the simple power 
law (12). 

Discussion. We have seen that the internal temperature Ti is determined 
by two temperatures, Tx and Ta, related to the non-degenerate and de- 
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log p [g cm-=] log p [g cm-3] 

Figure 1. Thermal profiles inside non-accreted (left panel) and accreted (right panel) 
envelopes of NS at two effective temperatures, logTe [K] = 5.9 and 6.5 (marked near the 
curves). Solid curves - numerical solution (Shibanov et al. 1998), dashed curves - analytic 
approximations (19) and (27). Straight lines marked "s", "d," and "t" give the values p 
and T at which the various temperature profiles cross the radiative surface [Eq. (21)], 
the onset of electron degeneracy [Eq. (22)], and the turning point [Eq. (24)], respectively. 
The Fe melting line is also shown. 

generate layers, respectively, and that the temperature growth occurs in 
the very surface layers of the star. This justifies the separation of the NS 
into a blanketing envelope and internal isothermal layers used in numerical 
simulations (GPE; PCY). We have also shown, in agreement with GPE, 
that the bulk of the envelope's thermal insulation arises in the relatively 
low density region around the turning point defined by Eq. (24). 

The dependence on Z and He that enters Eqs. (27) makes the Tj/Te 
ratio sensitively dependent of the chemical composition of the envelope. If 
we replace the iron envelope by an accreted envelope composed of hydrogen 
or helium, we get this ratio reduced by about a half order of magnitude, 
which corresponds to two orders of magnitude higher photon luminosity at 
a given internal temperature. Thus an envelope composed of light elements 
is much more transparent to heat, and this strongly affects cooling, as noted 
by PCY. Our analytic estimates present a tool for the fast estimation of 
the magnitude of such effects. 

In the next section we will see how these principal properties change as 
a result of a strong magnetic field permeating the NS envelope. 
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3. Effects of Magnetic Fields 

It is well known that the strong magnetic field can profoundly alter the 
physical properties of the NS outer layers ~ for reviews see Canuto &: Ven- 
tura (1977), Meszaros (1992), and Yakovlev & Kaminker (1994). We will 
be interested here in modifications introduced by the intense field B in the 
EOS, the Fermi temperature Tp, and the opacity K, all of which affect 
the heat transfer problem. Our primary focus will be on the properties of 
an electron gas. Free ions give only a minor contribution to the opacities, 
whereas their contribution to the EOS remains unaffected by the field in 
the non-degenerate regime and is negligible when electrons are strongly 
degenerate. 

Apart from the bulk properties of the electron plasma, magnetic fields 
also modify profoundly the properties of atoms which become very elon- 
gated and compact, having sharply increased binding energies (e.g., Canuto 
& Ventura 1977). This should strongly affect the emitted spectra from NS 
surfaces. Many works have been devoted in the past to the calculation of 
quantum-mechanical properties of atoms at rest in strong magnetic fields 
(e.g., Miller & Ncuhauser 1991). Some of them have been used to construct 
magnetic NS atmosphere models (Rajagopal et al. 1997). However, ther- 
mal motion of the atoms at realistic NS temperatures breaks down the 
axial symmetry and may completely alter atomic properties. It is there- 
fore necessary to have quantum-mechanical calculations for moving atoms, 
and their results should be included in models of EOS and opacities. Such 
models are available to date only for hydrogen (e.g., Potekhin et al. 1998, 
1999b, and references therein). For other species, this work still remains to 
be done. 

In the past, a lot of work has also been devoted to evaluating the con- 
ditions under which magnetic molecular chains may become stable in the 
surface layers; furthermore, in superstrong fields they may form a magnet- 
ically stabilized lattice (e.g., Ruderman 1971; Lai Sz Salpeter 1997). Such a 
phase transition is also expected to have observable consequences, but the 
field still remains largely unexplored. 

3.1. ELECTRON GAS IN MAGNETIC FIELD 

Motion of free electrons perpendicular to the magnetic field is quantized in 
Landau orbitals with a characteristic transverse scale equal to the magnetic 
length Om = {hc/eBy/'^ = X/Vb, where b = huJc/rrieC^ = Bu/M.li is the 
magnetic field strength expressed in relativistic units, lUc = eB/mc is the 
electron cyclotron frequency, and B12 = B/10^'^ G. The Landau energy 
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levels are 

/ 9 9 9A ^1'^ 

e = eniPz) = c [nigC + 2hucmen + p^j , (28) 

with n = 0, 1, 2, where the magnetic field B is assumed to be homoge- 
neous and directed along the z axis, and Pz is the longitudinal momentum. 
The ground Landau level n = is non-degenerate with respect to spin 
projection (s = —1, statistical weight go = 1) while the levels n > are 
doubly degenerate (s = ±1, gn = 2). The anomalous magnetic moment 
of the electron, g^ = 1.00116, causes splitting of the energy levels n > 1 
by 5e = {gc — 1)TilOc^ which, strictly speaking, removes the double spin- 
degeneracy. In typical NS envelopes, this splitting is negligible because 5e 
is smaller than either the thermal width ~ kT or the collisional width of 
the Landau levels. 

The electron's phase space is thus now a combination of an energy 
continuum in pz, corresponding to the motion along the field, and a dis- 
crete spectrum (the quantum number n) corresponding to the quantized 
transverse motion. This property will be refiected in most of the physical 
processes of our interest here. 

Let us denote by n* the highest Landau excitation populated at a given 
energy e. It equals an integer part of the combination pg(e) /(2mc/iajc), where 
Pn{f) = [{e/ c)"^ — {mec)"^ — 2mJ%UcnY^^ = \pz\- Taking into account that the 
number of quantum states of an electron with given s and n in volume V 
per longitudinal momentum interval Ap^ equals yAp^/(47r^a^7i), one can 
obtain the electron number density and pressure from first principles (e.g., 
Blandford &: Hernquist 1982). For strongly degenerate electrons, 

1 

= ^^E5nPn(eF), (29) 

n=0 



Pe = Po^Y.9n{l + 2bn) [xn^l + xl - ln(x„ + ^/l + xl)], (30) 

n=0 

where x„ = cpn{€p) / €n{0) , and Pq is the same as in Eq. (5). 

For a degenerate electron gas the thermodynamic quantities such as 
pressure, magnetization, and energy density exhibit quantum oscillations of 
the de Haas-van Alphen type whenever the dimensionless Fermi momentum 
reaches the characteristic values x = V^nb which signify the occupation of 
new Landau levels. In these oscillations the various quantities typically take 
values around their classical B = values, except in the limit of a strongly 
quantizing field (n* = 0) where one often finds substantial deviations (e.g., 
Yakovlev &: Kaminker 1994). The latter case takes place when the typical 
energies kT, kTp < ei — 1 - i.e., at T ^ and p < pB, where 

PB = m^UBlJie « 7045 /Xe g cm" , (31) 
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Tb = noJc/k-iY ~ 1.343 X 10^ (B12/7F) K, (32) 

and ub = l/('7r-^-v/2 Om) is the electron number density at which the Fermi 
energy reaches the first excited Landau level. This case is of special interest 
for the NS outer layers under consideration. 

Strongly quantizing field. When the electron's transverse motion is frozen 
in the ground state Landau level n* = 0, the phase-space is effectively 
one-dimensional. Then e = c (mgc)^ +p1i and Eq. (29) simplifies to 

pp/h = 27r^a^ne. (33) 

We therefore see that the dimensionless Fermi momentum, 

X = PF/meC = (2/3) xl/b « (0.6846/5) pe/^e, (34) 

is proportional to the density p, in sharp contrast to the non-magnetic 
Eq. (4). Henceforth wc denote ppo = i^cCXo the "classical" (non-magnetic) 
Fermi momentum at a given density, and reserve notation pp = itIqCX for 
the same quantity modified by the magnetic field. According to Eq. (31), 
the strongly quantizing regime in which Eq. (34) is valid requires x < V^b. 

The Fermi temperature Tp is again given by Eq. (3), but with the mod- 
ified X- Since x = i^/'^)^^^ (p/ Pb)'^^^Xo, Tp is strongly reduced at p <C pb- 
Conversely, at a given T < Tb, the degeneracy takes hold at much higher 
density than in the B = case. An initially degenerate electron gas at 
B = will thus become non-degenerate when a strong quantizing field is 
switched on. 

Let us now consider the EOS. Since n = 0, Eq. (30) simplifies consid- 
erably. Given Eq. (34), this expression again takes the form of a power-law 
of the density, P = -Po^X^'"*/(27r^7ad) oc p^'^'^ / B'^^~^ , with 7a,d = 3 or 2 in 
the non-relativistic and ultrarelativistic limits, respectively. 

3.2. MAGNETIC OPACITIES 

In the presence of a magnetic field, the conductivity k becomes a ten- 
sor, so that the heat fluxes along and across the field become different. 
Since the field varies over the NS surface, the heat transport becomes two- 
dimensional. Fortunately, since the crust thickness is relatively small, the 
one-dimensional equation (9) remains a good approximation, with n = 
K|| cos^ 6 + K± sin^ 9, where k\\ and k± are the conductivities along and 
across B and 9 is the angle between B and the normal to the surface. 



MAGNETIC NEUTRON STARS 



15 



Radiative opacities. In a magnetized plasma, two propagating polariza- 
tion normal modes are defined in the presence of an external field having 
widely different mean free paths each in the various photon-electron in- 
teractions (e.g., Meszaros 1992; Pavlov et al. 1995). Silant'ev & Yakovlev 
(1980) have calculated the Rosseland opacities for the cases when they are 
determined by the Thomson and frcc-frce processes. When Tjj ^ T, Kj- 
grows proportionally to {Tb/T)"^ - i.e., decreases as (T/Tb)'^- In partic- 
ular, at T <C Tb, the free-free opacities tabulated by Silant'ev &: Yakovlev 
(1980) tend to 

K,{B) « {23.2T/TB fK,{0) ~ 2.2 {Z/fiD pT-'-'' B^i cm^/g, (35) 

where K^{0) is given by Eq. (12). This estimate may be used if only 
Tq < Bi2, X 1) and the opacity is dominated by the free-free processes. 

Prom Eqs. (2), (20), and (35) one sees immediately that the radiative 
surface in the strong magnetic fields, B12 ^ Tq, is pushed to the higher 
densities, ps oa B. 

Caution is necessary however while using the scaling law (35). The 
strong magnetic field shifts the ionization equilibrium toward a lower degree 
of ionization, because of the increasing binding energies. Therefore, even if 
the plasma were fully ionized at some p and T in the absence of magnetic 
field, it may be only partially ionized at the same p and T when B is high. 
This increases the significance of the bound-bound and bound-free opacities 
and may result in a total radiative opacity considerably larger than that 
given by Eq. (35) (see, e.g., Potekhin et al. 1998). 

Electron conductivities. Unified expressions and fitting formulae for ther- 
mal and electrical electron conductivities in a fully ionized degenerate 
plasma with arbitrary magnetic field have been obtained recently (Potekhin 
1999). These conductivities undergo oscillations of the de Haas-van Alphen 
type at p > pb. At B ^ 10^^ G, the transport across the field is suppressed 
by orders of magnitude. This fact allows us to neglect k± totally, which 
will be a good approximation everywhere except at a narrow stripe near 
the magnetic equator, where 9 ~ 7r/2. 

This approximation, k ~ /ty cos^ 9, holds despite the arguments by Heyl 
& Hernquist (1998) that k± becomes non-negligible in the strongly quan- 
tizing regime at low density because the ratio evaluated for strongly 
degenerate electrons increases without bound as x ^ 0. However, the finite 
thermal width of the Fermi level (neglected by these authors) removes the 
divergency. At typical NS temperatures, the thermal averaging terminates 
the growth of k± and moderates the decrease of Ky, before they become 
comparable. In Fig. 2, we plot by solid lines an example of the thermal 
conductivities calculated according to Potekhin (1999). For comparison. 
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Figure 2. Longitudinal (||) and transverse (-L) thermal conductivities in the outer NS 
envelope composed of iron at T = 10* K and B = 10^'' G: comparison of accurate results 
(solid linos) with the classical approximation (dashed lines) and with the results without 
thermal averaging. The electrons are degenerate to the right of the vertical dot-dashed 
line, which corresponds to T = Tp. The dashed horizontal line in the non-degenerate 
region shows the conductivity given by Eq. (14) with Fz/A = 1. 



the dotted curves display the conductivities which would have been ob- 
tained without thermal averaging. As the electrons become non-degenerate 
at low densities, the solid curve for Ky is seen to level off. It tends to its 
non-degenerate value which is of the order of the non-magnetic value, 
Eq. (14), depending on density only logarithmically. 

The dashed lines show the "classical" approximation where the quan- 
tizing nature of the field is neglected. For Ky, this approximation is close 
to the non-magnetic one. We can see that at high enough densities beyond 
the first oscillation, the classical approximation is good enough. At lower 
densities, the quantizing nature of the field must be taken into account. At 
p < Pb, the neglect of thermal averaging is justified as long as the electrons 
are degenerate. Then the longitudinal conductivity may be written as 

^ k^TimjcHf X^^2_ Ko(l + xg)A , . 

'^11 127r^3^e4ne 2Q\\ 3 Qll ' ^ ' 

where kq is the non-magnetic conductivity given by Eq. (13), in which x 
must be replaced by xoj ^ is the non-magnetic Coulomb logarithm, and 
QII is a function of x defined by Eq. (A9) in Potekhin (1999). In the liquid 
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regime (far enough from the sohd phase boundary), the latter function 
reduces to the expression (Yakovlev 1984) 

Q" =r'-e«£;i(6, (37) 

where ^ = 2x^/6 + ^(ainQs)^! is an effective Coulomb-screening wave 
number, and Ei is the standard exponential integral. A simple order-of- 
magnitude estimate of «;|| in the degenerate Coulomb liquid can be obtained 
if we neglect qg and the second term in Eq. (37). In this way we obtain 

K|| ~ (4/3) kq A(1 + xo) X^/b ^ 5 X 10^^ Z'^ Tq erg/(cm s K). (38) 

As noted earlier, in the very strong fields considered here the onset of 
degeneracy is pushed into ever increasing densities as B increases. Therefore 
the turnover from radiative heat transfer to electron-conduction dominated 
transport may occur in the non-degenerate regime. In this case, Eq. (14) 
may be used to evaluate . 

3.3. CONSEQUENCES FOR THE HEAT TRANSPORT 

We have seen that in strong magnetic fields there are several different 
regimes regulating the EOS and opacities. For the construction of an ap- 
proximate analytic thermal profile it is sufficient to note that the non- 
magnetic expressions for the radiative opacity and longitudinal electron 
conductivity remain good approximations unless the field is strongly 
quantizing. Magnetic oscillations, which occur around the classical func- 
tions, will be smoothed out by integration while obtaining the thermal 
profile from Eq. (9); thus they are not too important. 

When the field is strongly quantizing, the opacities are modified appre- 
ciably. However, in the degenerate part of NS ocean, which is of our prime 
interest here, the analytic expressions for can be again approximated 
as a power law, Eqs. (38) and (14). The same is true with respect to the 
extreme quantizing limit of radiative opacity, as follows from Eq. (35). In 
the non-dcgcncratc regime, the magnetic field docs not affect the EOS. In 
this case, wc recover the solution (16) with new values of /3 = —1.5 and Kq, 
given by Eq. (35). Then 

n « 0.95 (5eff qf/' (p/Me)'/' B;^^/', (39) 

with the same q as in Eq. (19). Thus the temperature is reduced (its pro- 
file becomes less steep) with increasing B as long as the field is strongly 
quantizing {p < pb)- 

It is interesting to note that the value of the constant conductivity, 
Eq. (17) is independent of the magnetic field, while its numerical value is 
only slightly lowered as a result of the changed coefficient (3. 
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Sensitivity strip. As we have seen, the sensitivity strip is placed near the 
turning point from radiative transport to electron conduction, defined by 
K,. = K(.. In a strongly quantizing field, using Eqs. (35) and (38), we have 

p«i 250^e5eff°-'7^6-^|cosf?r°-^Si2 g cm'^^ (40) 

instead of Eq. (23) . With the temperature profile (39) , we now obtain 

_ 3.5 X 10' 6/31 10/31 K ~ 3 X 1(3^^ <f'-''- Bi2 3 

If, however, the electrons are non-degenerate along the turnover line, 
then K\\ is represented by Eq. (14) instead of (38). In this case, we obtain 
the turnover at 

pf^b2 sJk/Fz lie Te I cos 9\-^ Bu gcm'^ (42) 
Applying Eq. (39), we get 

2^X10^ /Ag\0-4 I.lxl03g04gl2 / A^°' _3 

Thus, in both cases (for degenerate and non-degenerate electrons) we 
have obtained similar dependences of the position of the turning point on 
the NS parameters. These dependences should be compared with Eq. (24). 
We see that Tt in both equations have similar values at ^ = 0, but in 
the magnetic field Tt increases with increasing 9. It is noteworthy that 
Tt is independent of B, while pt grows linearly with B in the strongly 
quantizing field. From Eq. (31) we can evaluate the condition at which pt 
lies in the region of strong magnetic quantization. Assuming that 9 is not 
close to 7r/2 and neglecting factors about unity, we see that pt < Pb for 
B12 ^ (2Tg^g/5i4)^^/^^, i.e., for the strong-field pulsars and magnetars. 

Onset of degeneracy. Let us estimate the point at which the electrons 
become degenerate. For simplicity, let us assume that the electrons are 
non-rclativistic. Taking into account Eqs. (3) and (34), we sec that the con- 
dition T = Tp in the strongly quantizing magnetic field is equivalent to p ~ 
608yLie\/?6-Bi2. Then from Eq. (39) we obtain p^ ~ 3700 {g^s q)^/"^ PeBu. 
Thus, similar to the non-magnetic case, the switch between the regimes of 
photon and electron heat conduction occurs not far from the onset of degen- 
eracy: Pd ^ Pt- Depending on 9, it can occur either in the non-degenerate 
domain (at 0) or in the degenerate domain (at 9 > 60°). 
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Figure 3. Temperature profiles through an iron envelope of a NS with M = 1.4 A/©, 
i? = 10 km, and with magnetic field B = 10^^ (left), lO" (middle), or 10^^ G (right 
panel). In every case, effective surface temperature was fixed to Te = 5 x 10® K (solid 
lines) or 2 X 10® K (dot-dashed lines). Different lines of each bunch correspond to different 
values of cosO: 1 (the lowest line of a bunch), 0.7, 0.4, 0.1, and (the highest line). 

Beyond the turning point. The integration beyond the turning point in 
order to obtain Ti can be done in the same way as in the non-magnetic 
case. However, since the turnover occurs in the strongly quantizing field, the 
integration path should be divided in two parts: below ps, where Eq. (38) 
should be used for the conductivity, and above p^, where one can use 
Eq. (25) with the right-hand side divided by cos^ 0. The result is similar to 
Eq. (27), but contains a profound dependence on the inclination angle: the 
thermal gradient grows rapidly as 6 increases toward 7r/2. 

This dependence is illustrated in Fig. 3, where we have shown several 
temperature profiles calculated numerically (Potekhin 2000). The curves 
start at the radiative surface, where T = Te, and end near the neutron drip 
point. The onset of the profound ^-dependence signals the turning point. 
According to our estimates above, this point is shifted to ever higher den- 
sities with increasing B. For the radiative boundary, the above-mentioned 
linear dependence ps oc S is confirmed. Another interesting effect is that T 
may continue to grow up to the bottom of the outer envelope, well beyond 
the "canonical" limiting density p = lO^'^ gcm~^. 

4. Conclusions 

In the preceding sections we have given a simplified analytic discussion of 
the neutron star outer envelope. It was possible in this simple picture to 
recover main features of more careful numerical calculations, such as the 
sensitivity of the temperature ratio T^/T^ to chemical constitution and the 
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general behavior of the (T, p) profile. 

Magnetic fields have a strong influence on the radiative and conduc- 
tive opacities in NS envelopes, as we have seen in the last section. It is 
remarkable though that in the region of the magnetic polar cap, in spite of 
the sharply reduced radiative opacities in the outer layers, the temperature 
ratios Tt/T^, and Ti/T^ are only mildly affected. The effect of the reduced 
opacity is partly counterbalanced by the increased mass of the radiative 
layer. Furthermore, when the magnetic field is inclined to the surface, the 
temperature gradient is increased because the electron conduction is effi- 
cient mainly along the field lines. Clearly more work is still necessary to 
fully analyze all these details. 

The possibility of magnctars, or extremely highly magnetized NS, in our 
Galaxy has been highlighted by recent observations of soft gamma repeaters 
(Kouveliotou 1998, 1999; Thompson 2000) and the special class of anoma- 
lous X-ray pulsars (e.g., Mcrcghetti 2000). Understanding the role of high 
magnetic fields in the observable properties of such objects has thus become 
an important theoretical task. Recent observations of compact galactic ob- 
jects further suggest rather small emission areas and high temperatures, 
which may be attributed to magnetars (e.g., Pavlov et al. 2000; Dar k. de 
Rujula 2000). The presence of an accreted hydrogen layer on the magne- 
tar's polar cap could allow for a rather sharp temperature contrast between 
the hot spot and the rest of the NS surface. This configuration has been 
invoked to interpret the observations in the case of the central compact 
object in the supernova remnant Cas A (Pavlov et al. 2000). 
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